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Abstract. We investigate the Collins effect in single spin asymmetries (SSAs) of the p^p ttX process, by 
taking into account the transverse momentum dependence of the microscopic sub-process cross sections, 
with the transverse momentum in the Collins function integrated over. We find that the asymmetries due 
to the Collins effect can only explain the available data at best qualitatively, by using our choices of quark 
distributions in the quark-diquark model and a pQCD-based analysis, together with several options of the 
Collins function. Our results indicate the necessity to take into account contributions from other effects 
such as the Sivers effect or twist-3 contributions. 
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Single spin asymmetries in hadronic reactions have become an active research topic recently since they may help us 
to uncover the transverse quark structure of the nucleon. Experimentally, there are data on meson electroproduction 
asymmetries in semi- inclusive deep inelastic scattering m21Elj whose most common explanation is to relate them to 
the transversity distribution of the quarks in the hadron ^BH], convoluted with a transverse momentum dependent 
fragmentation function |B1I7II51I^10| . i.e., the Collins function [Sllllj . The Collins function, which gives the distributions 
for a transversely polarized quark to fragment into an unpolarized hadron with non-zero transverse momentum, has 
aroused great interest recently since a chiral-odd structure function can be accessible together with another chiral-odd 
distribution / fragmentation function. 

Large single spin asymmetries (SSA) in p^p ttX have been also observed by the E704 Group at Fermilab |12j . 
Anselmino, Boglione and Murgia tried to reproduce the experimental data using a parametrization of the Collins 
fragmentation function, with several assumptions and in a generalized factorization scheme at the parton level. Notice 
that the Collins function obtained in Ref. is different from the original Collins parametrization 0. Other mech- 
anisms such as the Sivers effect |14II15II16II17) . and twist-3 contributions |18II19II2U) . have been also proposed. Very 
recently, Bourrely and Soffer [23 pointed out that the SSA observed several years ago at FNAL by the experiment 
E704 ^21 and the recent result observed at BNL-RHIC by STAR |22 are due to different phenomena. It is thus 
important to check whether a calculation with more detailed microscopic dynamics taken into account, and with also 
proper constraints on the Collins function, can indeed explain the single spin asymmetries in p^p ttX. We will show 
in the following that the Collins effect in single spin asymmetries (SSAs) of the p^p — * ttX process can only explain 
the available experimental data at best qualitatively, but we have not been able to explain the magnitude of the data 
at large a;^ by using several options of the Collins function and with our choices of the quark distributions in the 
quark-diquark model and a pQCD-based analysis. In our analysis, we work within the same simplified planar configu- 
ration of Ref. |13j . using several options of Collins function with a Gaussian- type transverse momentum dependence, 
together with our choices of quark distributions. Notice also that measurements of SSAs in weak interaction processes 
can distinguish between different QCD mechanisms 

We first focus on the formalism of description of single spin asymmetries at the parton level. For the inclusive 
production of a hadron C from the hadron A and hadron B collision process 

A + B^C + X, (1) 

the Mandelstam variables s, t, and u are written as 

s^{Pa+ Pb? = Ml + Ml + 2Pa ■ Pb , (2) 
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t=iPA- Pc? = Ml + - 2Pa ■ Pc, (3) 
M = (Pb - Pc? = Ml + 1/2 - 2Pb ■ Pc- (4) 

The experimental cross sections are usually expressed in terms of the experimental variables xp = IPi^j •Js and Pt 
at a given s, where Ph ~ and Pt arc the longitudinal and transverse momenta of the produced hadron C 

respectively, with energy Ec = \/ M'^ + P'^ = ^ A/^ + P| + = ^ il/^ + P^ + s/4 in the center of mass frame 
of the collision process, and s is the squared center of mass energy. The momentum and the energy of hadron A in 
the center of mass frame are 



(s - M\ ~ MlY - m\Mi 



PAL = \-PA\^\r- ^ f (5) 



Ea^^MI + PI^, (6) 
and the momentum and the energy of hadron B are 

Ps = -Pa, (7) 



Eb = ^MI + PI^. (8) 
The Mandelstam variables t and u can be expressed as 

t={PA- PcY = Ml + A/2 - 2EaEc + 2PalPl, (9) 
u^{Pb- Pc? = Ml + A/2 - 2EbEc - 2PalPl- (10) 

If the transverse momentum of the produced hadron relative to the fragmenting quark is negligible and the factorization 
theorem applies, the cross section can be written in terms of the subprocess at the parton level as 

dcr = -— = > / h/A[Xa)Jb/B[Xb)sd{s + t + u) -j (Xa, Xf,, z)Dc/c[Zj. (H) 

a,b,c,a 

However, there is no such a theorem if one considers transverse parton momenta as done sequently by taking the 
Collins effect into account. The Collins effect consists in taking into account the transverse momentum kc^ of the 
produced hadron C relative to the fragmenting parton c with momentum Pc, so we have 

Pc = 2Pc + kc±, (12) 

where z is the momentum fraction of the produced hadron relative to the fragmenting parton. The cross section is 
then expressed as 



dcT=^^._. 

(13) 



d'-'Pc 



EaAcM I d^k,^^^^-^fa^A{Xa)fb/B{xb)sS{§ + t + ^i) ^£±^1 (3;^ , , k,^)/?c/c(2, kci), 
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where D(7/c(-z, kc^) describes the fragmentation of hadron C with longitudinal momentum fraction z and transverse 
momentum kcj. relative to the fragmenting parton c. The denotations a, b, and c should include both the flavors and 
respective polarizations of the involved partons, if any or some of the incident hadrons A, and the produced hadron 
C are polarized. The cross section of the subprocess a + b ^ c + d 

do- = -T {Xa,Xt,Z,'kc±), (14) 

dt 

should be written in terms of the Mandelstam variables s, t, and u 

S = 2pa- Pb= XaXbS, (15) 
f=-2pa-Pc= ^t^'t(±fcc±), (16) 

u = -2pb-Pc='^u<Pu{±Ki_), (17) 
where the functions (I>t{±kci_) and ^ii(±fcc±) are given by [23 

^t{±k,i.) = <?(fce±) |l T 2fccx^|^ - [1 - <?(fcc±)] ^1 , (18) 

^u{±Ki_) = g(kc±) (i ± 2fcc±^^ + [1 - g(kc±)] ^1 , (19) 

1^ u(t + u) 2u ) 



with g{kc±) = -y/ 1 — k'^^ / , and ±fcc± referring respectively to the configuration in which kc_L points to the left 
or to the right of Pc. The kineniatical effect from the transverse momentum kc± is explicitly taken into account in 
<?t(±fcc_L) and <P-u{±kc±_), which become 1 for kc± = 0. The four variables Xa, Xb, and kc± are not independent, and 
by exploiting the S{s + t + it) function we get 

XaXbS + - [Xat<Pt{±kc±) + XbU<Pu{±kc±)] ^ 0, (20) 

z 

so we have 

^ ^ Xam{±kc±) + XbU<Pu{±kc±) ^21) 

After integrating the ^-function, we get the cross section 

(22) 

= Ea,hx,d/d^kc± dXa dXbfa/AiXa)fb/Bixb)i^ '^''''l^ {Xa,Xb, Z, k^^ ) | Dc/c{z, k^J,), 

where 



XbS + mi±kc±)' " S + U$u{±kc±) 

Because of the kc± dependence of x™*" and x™*" , the integration over kf,j^ should be performed after the integrations 
over Xa and Xb, and kc_L should meet the constraint kc± < Pc- Once we know the quark distributions fa/Ai^) smd 
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fb/six), the cross sections da ~ — {s,i,u), and the fragmentation functions Dc/^{z,'kc±), we can calculate the 
cross section (|^ explicitly. 

Now we apply the above formulas to the single spin asymmetries 

of ttX process, measured by the E704 Group at Fcrmilab After a calculation of the cross section in terms 

of the explicit polarization dependent ingredients, we find that the above asymmetries can be written as 

^(2; P ) - / <i'^'i^c±dxadxbSqa{xa)qb{xb)Aa{s, i, u)AD^{z, kcj_) ^^^^ 
/ d'^'kc±dxadxbqaixa)qb{xb)a{s,i,u)D{z,'kc±) 



where Sqa{xa) is the quark transversity distribution, qb{xb) is the usual quark distribution, da = da'^^ +da'^^, Ada 
do-Tf -do-T-'-, D{z,kc±) = [D^/c^{z,kc±) + D,,/ci{z,'kc±)]/2, and AD^{z,'kc±) = D^^/ci {z,kc±) - D^/ci{z,'kc±) is the 



Collins function. 

For the quark distributions, we use those in a quark-diquark model pHI and a pQCD based analysis |25) . explicitly 
taking into account the sea contributions based on the GRV parametrization of the parton distribution functions in 
Ref. |26j . The unpolarized and transversely polarized cross sections of the subproccss at the parton level, a{s, t, ii) and 
Aa(s, t, w), can be found in Refs. |27II28| . The detailed information on the quark transversity distributions in the quark- 
diquark model and the pQCD based analysis can be found in Ref. We need to point it out here that these transversity 
distributions do not come from model calculations, but from relations that connect the transversity distributions with 
parametrized unpolarized quark distributions, so that the calculation can be performed at the same scale at the 
experiment. It is important to use a same set of both unpolarized and polarized quark distributions, otherwise it 
is not reasonable to compare the denominator with the numerator. The calculated ratio would be unreasonable if 
the transversity distributions were taken from a model and then use parton distributions from a parametrization in 
order to perform the calculation. This aspect has been carefully considered in our calculation. In a strictly sense, the 
unpolarized quark distributions and transversity distributions evolve differently, so that we should use the relations 
at a specific initial scale such as « 2 GeV^, and then consider the evolution of the numerator and denominator 
separately. However, the effects of evolution arc presumably smaller than other uncertainties in the approach, such 
as the neglect of unfavored fragmentation, so we neglect the explicitly difference in evolution and take the energy 
scale the same in both the numerator and denominator. Besides, the quark-diquark model |24| has been successful in 
providing good descriptions of the quark hclicity distributions from experiments |29j , as well as data for nucleon form 
factors |3U| . so the extension of applying it to the transversity distributions is reasonable. 
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For the pion fragmentation functions, we take only the favored fragmentation into account, following Ref. and 
adopt the Kretzer-Leader-Christova parametrization [31] of D{z) 

D{z) = 0.689z-i O39(l _ z)i-24i. (26) 

The scale for D{z) is (Q^) = 2.5 GeV^. In our case s ~ 400 GeV^, but what matters for the evolutions is P^. Now 
Pl ~ XF\fs/'2,, which is about 5 GeV, and then P| is about 25 GeV^. We neglect the difference between different 
fragmentation functions in the evolution, as the effect due to the evolution, which is at most logarithmic, can be 
reasonably neglected in predicting the asymmetries when only ratios between different fragmentation functions are 
relevant. To take into account the kc^-dependence, we parameterize the fragmentation function D{z,\i.ci_) as [H] 

D{z, K^) - —D{z) exp(-i?2fc2^), (27) 

TT 

where R = . 2 with (k^i)^^^ = 0-44 GcV. The Collins function AD{z,kci_) is kc± odd, and it gives a null 

contribution in case kc± = 0, so we adopt the Collins parametrization for AD(z, kc±) as follows: 

AD{z,k,^) ^ 2Mc{kc^/z) 
2D{z,k,^) {AP^ + klJz^y 

where Mc is a typical hadronic scale around 0.3 1 GeV and we take its value as 0.7 GeV in our calculation. The 

above Collins parametrization fulfills the bound jl3j : 

\AD(z,k^i_)\ , , 

; ' /' < 1. 29 

2D{z,kci_) - ^ ^ 

In the above formulas, all quantities needed in the calculation are clearly given except that there are still large 

uncertainties concerning the Collins function AD{z, kc±). Therefore we consider several options for the Collins function: 

1. Option 1: The Collins parametrization Eq. H28|) . 

2. Option 2: The Collins parametrization with an additional z-dependent factor 

^Diz,k,^_) _ _2MciK±M_ 

2D{z,k,^) {Ml + klJz^f^'^' ^^"^ 
where P{z) = cz"(l — zY with a = 15, /3 = 1.5, and c ~ 306.18. This option is chosen to illustrate a case which is 
of similar shape and not far from that of the Anselmino-Boglione-Murgia (ABM) parametrization (Option 3), but 
which can reproduce the asymmetries with magnitude around twice that of Option 3. 

3. Option 3: The Anselmino-Boglione-Murgia (ABM) parametrization (there is some trivial difference in details) ^31 

AD{z,k,^) _ { 0.007409z(l - z)-^-^ z < 0.97742; 

[ 1 .> 0.97742, 

1 /'^ 

at an average value of (k^j^) " = fcj?j^. 
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4. Option 4: The upper limit bound 

2D(z,fc,^)-'- ^^^^ 

The z-dependence of the ratio '^^(^^'1^"^^ at a given kc±_ = = Ale for the four options can be found in Fig. ^ 
from where we can find tlrat the shapes of the Collins function are quite different for the four options. 



1 .2 




Fig. 1. The z-dependence the ratio '^^^^^'^ "^^^ at a given value of kc± = k'^j_ = A'lc for the four options of the Collins 
function. The solid, dashed, dash-dotted, and dotted curves are corresponding to the results for the 4 options of Collins function 
respectively. 



Thus we have all of the ingredients for the calculation of the single spin asymmetries A{xp, Pt) in (|25|l . The results 
shown in Fig. |5] are calculated at Pt = 1.5 GeV with the quark transversity distributions from both the quark-diquark 
model and the pQCD based analysis. We notice that both models give large A{xf,Pt) at small xp, with magnitude 
compatible with or larger than the data, but at large x p the magnitude is below the data. The quark-diquark model can 
reproduce the trend of A^xp.Pp) for tt^ qualitatively, but the pQCD based analysis produces a result with opposite 
sign. The reason for the discrepancy is that the valence o?-quark transversity distribution is positive at large x in the 
pQCD based analysis j9j. The situation for the quark-diquark model is better, but there is still a discrepancy with the 
magnitude of the data. In order to reduce the discrepancy between the calculated result and the data, especially for the 
TT^ data, we introduce an additional case (which is not a realistic case, because it breaks the Soffer's inequality |32| . so 
this is just as an illustration) in the quark-diquark model with the valence c?-quark transversity distributions enhanced 
by a factor of 3, i.e., 

8dv[x) = ~dv{x)Wv{x), (33) 

where (x) is the unpolarized quark distribution and {x) is a Melosh-Wigner rotation factor to reflect the relativistic 
effect from quark transverse motions |33) . We find that the magnitude of the calculated results is more compatible 
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0.5 1 0.5 1 0.5 1 

Xp Xp Xp 



Fig. 2. The single spin asymmetries of p^p ivX process. The calculated results in the quark-diquark model: (a), (b), and (c), 
in the pQCD based analysis: (d), (e), and (f), and in the quark-diquark model with more negatively polarized valence d-quark: 
(g), (h), and (i). The solid, dashed, dash-dotted, and dotted curves are corresponding to the results for the 4 options of Collins 
function respectively. The experimental data are given by the E704 Group |12| 



with the data than both the quark-diquark model and the pQCD based analysis, as shown in Fig. [3 This might 
suggest that the d-quark transversity distributions is more negatively polarized than predicted in the quark-diquark 
model, as supported by sum-rule based arguments j34| . if one only considers the Collins effect. From the results of 
the upper limit bound of the Collins function, we notice that any improvement in the parametrization of the Collins 
function cannot improve the fit to the data at large xp. Inclusion of unfavored fragmentation will lead the calculated 
results to go in the opposite direction from the data. The ABM parametrization reproduces the shape of the data, 
but underestimates its magnitude by a factor of around 4-6. We have neglected the gluon transverse polarization in 
our calculation, as the gluon transversity in a spin 1/2 hadron is strictly zero due to helicity conservation. Thus the 
description of the data is not obtained in our calculation by taking into account the detailed microscopic cross sections 
with the transverse momentum in the Collins function integrated over, together with our choices of quark distributions 
and Collins function. 
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In summary, we checked the Collins effect in single spin asymmetries (SSAs) of the p^p ttX process by taking 
into account the transverse momentum dependence of the microscopic sub-process cross sections at the parton level, 
with the transverse momentum in the Collins function integrated over. We introduced several options for the Collins 
effect, and found that the single spin asymmetries of p^p ttX process due to Collins effect can only explain the 
available data at best qualitatively in some specific situations, by using the quark distributions in the quark-diquark 
model and in a pQCD-bascd analysis. The results suggest the necessity of taking into account contributions from other 
effects such as Sivcrs effect and twist-3 contributions. It might be also possible that some unexpected novel behaviors 
of the quark distributions and Collins function need to be introduced. With our present calculation results on the 
asymmetries and the calculation done in Ref. for the cross sections themselves, we have to conclude that one needs 
to introduce other mechanisms in order to understand the single-spin asymmetries for tt inclusive production in pp 
collisions. A similar conclusion with a more detailed and complete analysis has been also drawn in a recent work |35j . 
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